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The result from the SNla projects suggest that the dark energy can be represented by a fluid 
with w < — 1. However, it is commonly argued that a fluid with \w\ > 1 contradicts causality. 

Here, we will show that a fluid with \w\ > 1 does not contradict causality if w is not constant. 
Scalar field are the most promising candidates for describing the dark energy and they do not have 
a constant equation of state parameter w. Scalar potentials may lead to regions where w is larger 
or smaller than one and even regions where the group velocity dp/dp diverges. 

We study the evolution of scalar field perturbations and we show that the "sound speed" is always 
smaller than the speed of light independently of the value of w = p/p or dp/dp. In general, it is 
'^f-' ■ neither the phase velocity nor the group velocity that gives the "sound speed". In the analysis 

we include the special cases of \w\ > 1, tachyonic, phantom, k-essence scalar fields and curvature 
corrections. Our results show that scalar fields do not contradict causality as long as there is 
■ dispersion. 

The existence of a dark energy, an energy density p with negative pressure p and equation of state w = p/p, has 
been determined by the SNla [2| and the CMBR observations [lj. These observations show that we are living in a flat 
universe with a matter contribution today Q m ~ 0.3 and a dark energy Qde — 0.7 with Wde < —2/3. The cosmological 
results do not rule out an energy density with w < — 1 . 

Perhaps the most economic solution to the dark energy is that of a scalar field. Not only because they are widely 
predicted by particle physics but also because their dynamics may lead in a natural way to an accelerated universe. 
These scalar fields are called quintessence and they are homogeneous in space and have only gravitational interaction 
, with all other fields. Other interesting possibilities have been consider such as phantom fields Q, k-essence Q or 
£NJ tachyonic scalar fields and curvature corrections |icj-[l^. 

t-H ' The possibility of having a fluid with w < —1 is problematic. There are objections against such a fluid from 
vacuum instability, causality and sound speed among other considerations. For a fluid with a constant equation of 
state parameter w it cannot exceed the value of \w\ > 1. One of the many problems appears as the propagation of 
sound, since for w constant it gives the sound speed and clearly w must be smaller than the speed of light c = 1 ||. 
Another related problem is the vacuum instability for fluids with a constant w < -1 0. Here, we are interested in 
studying the problem of causality and sound speed for fluids with \w\ > 1 or \dp/dp\ > 1. As we will see in the next 
section a scalar field with \w\ > 1 and/or \dp/dp\ > 1 can be easily obtained if the field has negative potential [2(j 
| (e.g. a tachyon, m 2 < 0) or for a phantom field with negative kinetic energy and positive potential 7]. Concerning 
the evolution of the perturbations of the scalar field we will show that in the absence of gravity a tachyonic field is 
| equivalent to a phantom field. 

We will show that the information speed for a scalar field does not contradict special relativity (see ^(|) even in 
the cases where \w\ > 1 and \dp/dp\ > 1. However, in these cases the amplitude of the perturbations grows rapidly 
and therefore the first order approximation used to study the sound speed fails. This signals that the space (universe) 
is no longer homogenous and non-perturbative mathematical methods must be used to analyze the behavior of the 
universe. Still, we want to emphasis that scalar fields with \w\ > 1 and \dp/dp\ > 1, tachyons, phantom or k-essence 
models do not violate causality and the sound speed is smaller or equal to the speed of light. 
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I. MOTIVATION 



It is usually stated that the "speed of sound" for a fluid is given either by = p/p if it is constant or by dp /dp 
if Wfj, is time dependent, where p is the pressure and p Q the energy density. In the case that w = p/p is constant then 
the two quantities dp /dp — w coincide. The requirement that the propagation of information should not exceed the 
speed of light implies that \w^,\ < 1 and/or \dp /dp \ < 1. However, as we will show, neither of these two conditions 
are necessarily satisfied for a scalar field. 
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Considering an homogenous scalar field, it is well known that its evolution leads naturally to a non-constant uu^, 
so wj, should not be interpreted as the " speed of sound" . The equation of state parameter for a scalar field is 

V4,o \\4>l~V{<j) ) 
w<h = - L - — —r-. 1) 

P0 O \\4>l + v(cp ) 

where the o denotes that we are considering an homogenous scalar with <j> (t) and 4> = d<p/dt and we have included 
the parameter A = ±1 to take into account a canonical scalar field A = 1 and a phantom field A = — 1. It is easy to 
see from the above equation that for any scalar field with positive kinetic term Ek = <p 2 /2 > that the magnitude 
and sign w<j, depends on the value of the potential V((f>) and eq.Q) gives 

Ko| < 1 V(<f> ) > (2) 

w^o > 1 ^ E k > -V > 
w^ < -1 & E k + V < 0. (3) 

As we see from eqs.© we can obtain w^o < — 1 for a scalar field with a negative potential (e.g. for a tachyon field 
with V — m 2 <j) 2 , m 2 < 0). We can also have w^, < —1 for a phantom field with negative kinetic energy (Ek < 0) 
with positive potential (V(4> ) > 0). It was shown in [2(| that this kind of behavior is reached naturally for scalar 
field with a negative minimum in the presence of a barotropic fluid that can be, for example, matter or radiation. 
Therefore, we do not need exotic fluids to have regions with \w,p\ > 1. 

Now, let us analyze the quantity dp^/dp^ for a scalar field with arbitrary potential. The homogenous part of the 
scalar field depends only on time and using eq.Q we have the adiabatic quantity 

d i± = t± = l 2V ' 1 1 2 YL (4) 

dp4, f>4, \4> + V 3A#c/> 

where we have used 

P<p = ^ - V') p = 0(A0 + V) (5) 

and the equation of motion of the scalar field 

X(4> + 3H<p) + V' = (6) 

for the second equality of eq.(@J with V' = dV/d<fi. Notice that in general eq.Q) has not only regions where {dp,}, / dp<j,\ 
is larger than one but regions where it diverges. This will happen at the turning points of the field oscillations 
around the minimum of the potential when <j) = 0, if the point is reached at a non-extremum of the potential V, (i.e. 
\(f> = —V' 7^ 0), e.g. for any potential with a minimum at finite <j> as for a massive scalar field with V = m 2 4> 2 . 

Clearly, we see that under normal conditions both (homogenous) quantities w^, dp^/dp^, can be larger than the one 
(the speed of light). Does this imply that the sound speed of scalar fields travel faster then the speed of light? The 
answer is no, as we will show later, and the reason is that neither w§ nor dp^jdp^ gives the correct interpretation of 
the sound speed. 



II. SUMMARY OF DIFFERENT VELOCITIES 



In order to clarify the notion of causality for fluid or a scalar field it will be useful to define the different velocities 
involved in determining the transportation of a signal in a medium. 

For simplicity in presentation purposes we will consider only one space dimension. The wave equation for a fluid is 
given by 

v 2 

p- -r- dlp + d p + bp = (7) 

&■ 

where c is the speed of light and the parameters v, b, d may depend, in general, on the coordinates t, x and they define 
the type of medium the fluid is in. If we take a Fourier transformation 

/oo 
dk P k(i)e lkx (8) 
-OO 
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cq.ljjj takes the form 



Pk + d Pk + (—k 2 + b) p k = 0. (9) 



The parameter v defines the phase velocity while the friction term is given by d. The simplest wave equation is 
given by d = b = and v constant. In this case all monochromatic wave functions will have the same speed and 
v can be interpreted as the sound speed. For d = and (v 2 k 2 /c 2 + b)/k 2 is fc-dependent, the medium is said to 
be dispersive and the monochromatic wave functions will travel with different phase velocities. In this case it is the 
group velocity that determines the speed of sound. If we allow d to be different than zero then we have a dissipative 
medium and the wave amplitude will suffer a damping or growth depending on the sign of d and neither the phase 
nor the group velocity would give the " sound speed" . 

For d = b = and v — cte the solution to eq. © involves a sum of monochromatic waves 

p(t, x) = f U k e - t(kx -^ + B k e Kkx-^t)\ dk ( 1Q ) 
J —oc 

with tu, k constants and the amplitudes A k , B k also constant. From eq.10, with b = d = 0, the frequency zu (it should 
not be confused with w the equation of state parameter defined in the previous section) is a function of k 

vo 2 = t,k 2 . (11) 

The wave number is given by k = 2ir/\ and gives the inverse of the wave length. From now on we will set the speed 
of light back to one (c = 1). The speed for a monochromatic wave, e.g. the speed of the maximum of the wave, is 
given by the constrain kx — zut = cte from which we have the phase velocity 

dx zu 

v ^lu=J = v (12) 

if k,zu are independent of x, t. A velocity smaller than c requires v p — zu/k < 1. Notice that in this case v p is the 
same for all monochromatic wave functions since w/k is fc-independent. 

If we are in a dispersive medium 6^0 and the wave frequency depends on the value of k. In this case one needs 
to consider the evolution of wave packets since their shape will be time dependent. A wave packet occupies a limited 
region in space. It is common to assume that the monochromatic frequencies that form the packet are concentrated 
around a central value to with corresponding wave number k. For a wave function with a slightly different wave number 
k' = k + Ak the corresponding frequency is zu(k + Ak) and for Ak small we can write zu(k + Ak) ~ va(k) + (dzu/dk)Ak. 
The solution to the wave equation for a wave packet is of the form 

p(t, x)=A e^ kx -^f((Akx - tAvu)) (13) 

and the amplitude of the wave is modulated by the function /. The group velocity is defined by requiring the argument 
of / to be constant, i.e. by the solution of Akx — tAzv = cte, giving the group velocity 

Azu dm , 
V °=Ak = lk (U) 

if dza/dk is t, x independent. The group velocity v g is valid for any dependence of tu(k) but only as long as Ak is small. 
If the medium is non dispersive vu is constant then the group velocity and the phase velocity coincide v g — v p = w/k. 

In general the form of the packet will not remain constant and it will be smoothed out during its propagation and 
the condition of having a small Afc will no be maintained. 

Lastly let as consider a dissipative medium in the simple case of constant b, d, v. In this case the solution to eq.10, 
with the ansatz p k — Ae mt + Be~ ltnt give the equation, 



v 2 

zd 2 +idw+ ( — k 2 + b) = (15) 



with 



id±JA(Kk 2 + b)-d 2 
*>= (16) 
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FIG. 1: Propagation speeds near an absorption band. Although the phase velocity becomes larger than light speed and the 
group velocity becomes infinite and even negative, the signal velocity is always smaller than c. 

and we see that w takes imaginary values if d ^ or 4(k 2 v 2 jc 2 + b) — d 2 < 0. This will give damping or a growing 
wave function solution. If 4(k 2 v 2 /c 2 + b) — d 2 > the damping term will be the same for all monochromatic modes. 

The short explanation we have given above also shows conditions for the validity of the approximation: the packet 
must be quasi-monochromatic and the frequency a slowly varying real function of k. If either of these approximations 
breaks down, the group velocity loses its physical meaning. In the general case the "speed of sound" is neither 
the group nor the phase velocity and one must use the more complex and subtle notion of signal velocity [25l l2tS| . 
The signal velocity is defined as the velocity with which a given standard amplitude of the wave packet moves, for 
instance half that of the maximum amplitude. The standard Sommerfeld-Brillouin definition amounts to this. The 
condition of retarded wave guarantees then that the signal velocity will be always smaller than c, even if the wave 
packet deforms very much and the group velocity becomes v g > c. That this may happen is shown in the classical 
Sommerfeld-Brillouin diagram [22| for the three velocities in the neighborhood of an absorption line (Fig. In the 
case of a scalar field the group velocity diverges when = with V ^ 0. This will happen in most cases and even 
for a potential of the form V = m 2 (f> 2 around the turning points of the scalar field. 

III. PERFECT FLUID AND SCALAR FIELD 
A. Perfect Fluid 

A perfect fluid has an energy momentum tensor T^ v = (p + — pg^ where p is the energy density and p 

the pressure of the fluid, g^ v the metric tensor while the 4- velocity. In the rest frame u° — l,u l = and T^ v 
gives the usual energy momentum tensor T^ v — diag[p, —p, —p, —p]. For an arbitrary velocity v the four velocity is 
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= u°(l,v % ) with u° = 1/vl — v 2 . The equation of motion are given by the conservation equation 9 M T Ml/ = and 
to solved them one needs to add an equation of state p = wp. Matter and radiation have a constant w equal to zero, 
1 /3 respectively. 

In order to study the propagation of sound it is common to take perturbations around a centra! value p Q1 p Q 

p(t,x) =p {t) +p 1 (t,x), p(t,x) = p (t) + pi(t,x) (17) 

where the central values do not depend on the space coordinates while the perturbations p\, pi depend on t and x. 

Keeping first order terms in the pressure and energy density perturbations one gets by solving the equation of 
motion, with a Minkowski metric, the equations 

% = -(p + p)V, g = — ^Vp, (18) 
ot at p + p 

Using the relationship p = wp (or equivalently p\ = wp\) one has 

Vpi = wVpi (19) 

which is valid if w is not a function of x. The restriction on w — P\j P\ = dp/ dp to not depend on the spatial variations 
follows from the hypothesis of adiabatic perturbations. Eliminating v from eqs.lJTS} and using (|T§|| we get the usual 
wave equation 

pi - w 2 V 2 pi = 0. (20) 

Comparing ea. ((20() with eq.!© we see that the phase velocity is given hy v p = \w\ = v and it is the same for all 
monochromatic wave functions of the fluid perturbations pi since b = d = 0. Ea. 1)20(1 is valid for a static (Minkowski) 
space in the absence of gravitational interactions. 

If we consider a flat FRW space, with metric ds 2 = dt 2 — a(t) 2 dx 2 , then the perturbation equation for radiation 
with S r = Sp r / p r is given (in the transverse gauge) by 

Sr + ^k 2 6 r = -S c (21) 

where S c = 5p c / p c is the perturbation of (cold) matter and the r.h.s. term in ea. 1(21(1 is a source term for S r . Comparing 
eq. (|2"T|l and eq. we see that the phase velocity is given by v 2 — 1/3 and in the absence of the source term, the group 
and the phase velocity would coincide. Since there is no dissipative term (d — 0) the sound velocity for radiation 
would be given be v p = v g — v p = 1/ v3- 

Notice that even in an expanding universe the perfect fluid (in this case for radiation) has no dispersion. 



B. Scalar Field 



A scalar field is defined by the action 



S = dx 



(22) 



where g = det[g^] and we have introduced the parameter A = ±1 in order to incorporate phantom fields in the 
analysis (+1 for a canonical scalar fields and —1 for phantom fields). From this eq. we can extract the energy 
momentum tensor T^ v ((f)) = Xd^d^cj) — g^ v \^\g a ^d a <\>dj3(\) — V (</>)] and the equation of motion is 



A(d"d M <£ + + V = 0. 



(23) 



where we have set [d l i{g tlv \J—g)/\/—g\d li <t) = 3i/0, valid in a flat FRW metric, with = a {t) 3 i H = a/a the 

Hubble parameter and a(t) the scale factor. The T° — p component defines the energy density while T\ = —p<f, the 
pressure, for simplicity and presentation purposes we will only consider perturbations in one direction i — x, giving 



- - -aw) + v(<p) 



p<t> = k\^ 2 



(24) 
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and the equation of state parameter is 

P4> A(±0 2 - ±5^) + ^) 

In cosmology it is usual to assume the scalar field to be homogenous in space and therefore (f>(t, x) = 4> (t) is only a 
function of time and the space derivatives vanish, <9j</> = 0. In such a case ea. (|25)) reduces to eq.Q). 

However, contrary to a perfect fluid, in the case of a scalar field the equation of state cannot be imposed and it is 
dynamically determined by the evolution of (j> given by ea. (|23[l . It is easy to see that w in ea. l|25[l does depend on the 
spatial coordinates. Taking a perturbations in the x-direction around the average solution cj) one has 

<j>(t,x) = <f> (t)+Lp(t,x) (26) 

and to first order in ip we have V(<f>{t,x)) ~ V {cj> {t)) + V> + 0{ip 2 ) with V' Q = dV/dcf)]^ and (</> 2 - d l <j>d l <j))/2 = 
4> 2 /2 + (p ip + 0(ip 2 ). Ea. lt231) up to first order in ip gives 



2A VJ<t> 2 <p(t,x)-2V (t> <p(t,x) 

Wtj, = W^o H j 



t% 2 ) (27) 

where u^ is given by eq.Q) and p<p — X<p 2 /2 — V ((j) ). The perturbed equation l|23[) up to first order in ip is 

X(p - V'V + 3Hip) + V{<t> )"<p = - X -4> h (28) 

where the r.h.s in eq.JSHJ) is due to the variation of the metric tensor = rj^ + and V' 2 = a _2 V 2 and 
V" = d 2 V/d(/) 2 . In terms of a Fourier transformation <p(t,x) — (2n)~ 3 J dk 3 pke lk ' x ^° we can write ea. (!28ll as 

.. , „„. , ,fc 2 , V(Ut)Y\ 1 ; ; 

ipk + oHtpk + H : )<pk = --^TVkoh (29) 

a -2 A 2A 

with k 2 — k ■ k. Ea. (|29|) gives the wave equation for the scalar field perturbations and we can see that the medium is 
dissipative (H ^ 0) and dispersive ((k 2 /a 2 + V" /X)/k 2 is k dependent). In conformal time adr = dt ea. (|29|l becomes 

tf? + 2Hp% + (fc 2 + a2y{4> ° {t)) " ) Vk = ~4iLh* (30) 

where * means derivative w.r.t. dr and H = a* /a. Comparing eci. (|3U[) with eq.@ we see that v p = 1, i.e. the "phase 
velocity" is the speed of light, b = a 2 V{4> {t))" / A and d = 2H are non zero and time dependent. Therefore, the wave 
function will have a non constant amplitude and the velocity of each monochromatic wave will differ thus giving a 
group velocity different than the phase velocity. We also remark that the left hand side of eq. I|3U|) is the same for a 
canonical scalar field with negative mass square (i.e. a tachyon field with V"/X = V" < 0) to a phantom field with 
positive m 2 (i.e. V"/X = —V" < 0). On the other hand, the gravitational contribution given by the r.h.s of ea. H30|) 
has opposite sign for a canonical field and a phantom field since A changes sign. 

C. K-essence 

Let us now consider the following (k-essence) action for a scalar field Q 

S = J dx A ^g~K{cj))F{X) (31) 

with K((f>) an arbitrary function of the scalar field 4> and F(X) a function of the kinetic term X = <9 M 0<9 M 0/2. The 
equation of motion of ea. l|31fl . in a flat FRW metric, is given by 

(a(t) 3 F'(X)K(^d^) - a(t) 3 F{X)K'(cp) = (32) 

where we have taken yj—g = a 3 (t) and a prime denotes derivative w.r.t. the argument and we used SX/Sd^cj) = d^cj). 
Eq.(22J) is then 

d^<j> + + ^d^Xd^ + 3Hj> - = (33) 
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The Lagrangian in eq. I|31(l can be reduced to a canonical kinetic term scalar Lagrangian for F = X, K = 1 and eq. 1331) 
reduces to d^d^cj) + 3H<fi = 0, i.e. a scalar field without potential. For an homogenous scalar field 4> = 4> (t) we have 
d^X — <j) (j) and ea. (|3~3")l becomes 



d^fa + 3H4 ——r- + -— 1° . = (34) 




In order to analyze the evolution of the perturbations we will expand ea. (|33|l to first order in <j>(t, x) = <fi (t) + ip(t, x). 
Let us defined the functions D(<f>) = K'{<t>)/K{4>), E(X) = F"(X)/F'(X) and S(X) = F(X)/F'(X) which appear in 
eq.®- We expand the functions D, E, S to first order in <p(t, x), i.e. D = D + D'p 7 E = E + E'SX, S = S + S'5X 
with SX — 4> (p where the prime denotes derivative w.r.t. the argument. The expansion of X and d^X are given by 

x = x + sx = x + 4> a ip + o{ip 2 ) 

d^X = d^X. + ipd^o + ^d^ + Oiip 2 ) (35) 

with X = <frl/2, X = <j>o4>o, d a X Q = d a (d fl 4>od tl (f) /2) = (f> (j) 8oa and d^Xd^cf) = (j) X + 2tj) (f) ip + § 2 tp. To first 
order in <p(t,x) eq.|J53J| becomes 

<p-BV n (p + WH<p + A(p + m 2 eff <p = --<j) h (36) 

where 

B = 1/(1 + E <j%) = F'/(F' + F"tj} 2 ), A = B{2D ^ + E'4?J - S'D ^ + 2Ej <j> ) (37) 



and the effective mass is given by 



il ff =B(j>t-S )D> = ^ f°jf . (38) 



Eo^l 



Making a fourier transformation of oa . ( f36|> . as in ca. (|29|l . we get 
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(p k + (3BH + A)ip k + (-2 B + m 2 eff )<p k = -|0 o /i (39) 

The r.h.s. in ca. (|36|l is due to the variation of the metric as in ea. (|28|l and the factor of B in eas. l3(il) and (|38|l arises 
because d^Xd^cj) involves a second order time derivative of cp. Comparing equation (|39() to ea. (|29ll we see that they 
are similar, ea. l|39|l has a second time derivative term, a friction term and an effective mass fn\ff- The friction term, 
the effective mass and the phase velocity v 2 = B, as seen from eq.©, are model and time dependent. The phase 
velocity will, in general, be different than the speed of light contrary to a canonical scalar field. The effective mass 
square can in principle be positive or negative and the effect of the friction term is to damp or enhance the amplitude 
of the fluctuations, depending on its sign. 

Let us now consider the particular case presented in [8|. In this special case K{4>) = l/cf) 2 and F{X) = g(y)/y with 
y defined as y = l/Vx and the pressure, energy density, equation of state and group velocity are given by i, 

p = — — ( 4 °) 



TV 

g(y) _ g{y) - yg'iv) 



yg'{y) ' y 2 g"{y) 



(41) 



As seen from eas. (|41f) the equation of state parameter (i.e. the phase velocity) diverges when g'(y) = while the 
group velocity diverges at g"(y) = 0. In this class of models we have D(<p) = K 1 '(<p) / 'K(<p) — —2/<fi, D' = dD/dcj) = 
2/ti, E(X) = F"(X)/F'(X) = y 2 (y(g' + yg") - g)/(2(g - yg')) and S(X) = F(X)/F'(X) = 2g/(y 2 (g - yg')) and 
the effective mass of ea. (l38H is given by 

m 2 4 [<j>ly\g-yg')-2g\ 
eff y 2 €[2(g-yg') + (y(g + yg')~g)y 2 <P 2 ]' 

The coefficients A,B of ea. H37l) can be easily obtained. However, even in this particular case the effective mass, A 
and B have not a simple form and are still model and time dependent but the main conclusions remain valid. The 
sound speed for k-essence differs from a canonical scalar field, its phase velocity is no longer equal to the speed of 
light but and its sound velocity is smaller than the speed of light. 
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D. Curvature Correction 

In the case where we consider corrections to the Einstein general relativity we have a lagrangian [T3L ITij 

S = [ dx 4 V=g(f(R)+L rn ) (43) 



with f(R) an arbitrary function of the Ricci scalar R and L m (<j)) the matter lagrangian. The equation of motion of 
ea. (|43|l in a flat FRW metric, is given by 

f'(R)R aP - \f{R)g a(i = f(R)^ {g^gpu - g a pg M u) + f% (44) 



which can be written in a suggestive form [l3| 



1 



where 



R aP - ^Rg a0 = + T™p (45) 
T% = j^ ) \^[f{R)-Rf'(R)] + f\R)^(g a ,gp v -g aP g^ (46) 

T ^ " jwf" (47) 

In some interesting cases, such as 

f(R) = R§, f(R) = R-£ (48) 

these fourth order equations can be transformed to the Einstein frame, where the gra vitational equations take the 
Einstein form, coupled to a scalar field. This is done through a conformal mapping [13l Il4| 

g a p = F(R)g a0 , = G(R) (49) 

with suitable chosen functions F(R) and G(R). The resulting equations are 

R aP -\R~g a p = T+t + T^ (50) 
V M V0-V"(0) = a (51) 

for a suitable source term a |l3lll4| . Thus, in many interesting cases, curvature corrections model of dark energy can 
be mapped on the scalar field quintessence paradigm. Therefore, the small disturbances equations will have the form 

It should be noted that if the action i|43|) is treated with a Palatini variation 0, ^| the resulting field equations are 
of the second order type. The resulting field equations have an effective cosmological constant which, for a suitable 
choice of parameters, can be fitted to reproduce the present deceleration parameter. Since these equations have 
the Einstein form, small perturbations will behave much the same way as in a cosmological model with a varying 
cosmological constant. In particular, adiabatic scalar modes will obey a scalar ave e quat ion, formally similar to l|29|l . 



quati 

i la. 



However, for lagrangians singular in R — the newtonian limit cannot be recovered 

IV. SOLUTION 

We shall study now the propagation speeds of a scalar field. For the sake of simplicity, let us consider Eq. (I23|) 
in Minkowski space-time, with \J—g = 1 in a 1+1 dimensional space-time. Let us consider the Lagrangian density 
L = ^{(d c t(t)) 2 — (d z (j)) 2 ] — V((j>), where the constant A allow us to include phantom fields in the following discussion. 
The equation of motion as given by ea. (|23|l and the perturbed equation of motion for (p(z,t), ea. (|28|l (4>(z,t) — 
<fro(t) + ip(z, t)), in the absence of the source term, is 

1 d 2 w d 2 w\ 9 , , 

-f +mV = 0. (52) 



c 2 dt 2 dz 2 
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where we have taken V(<p )" — m2 ■ ln terms of the Fourier transformation we have 



A 



~2<Pk+ (Afc 2 + m 2 ) ip k = 0. 



(53) 



We have not incorporated a friction term because it only enhances or damps the amplitude of the fluctuations but 
does not interfere with the sound speed. Eqs.JH2l or tBi} encompasses several particular cases commonly found in 
the cosmological literature: a standard scalar field A > 0, m 2 > with 1 > w^ > — 1, a negative potential field 
A > 0, to 2 < and a phantom field A < 0, m 2 > 0. The last two examples have regions with \w^ \ > 1 and/or 
\dp /dp \ > f . The second case corresponds to <j) near a maximum of the potential, which means the existence of an 
effective tachyon. It would seem that the sound velocity becomes grater than c in this case, but this is not so: the 
fluid model breaks down when the propagation of perturbations in the system is considered, even in the "normal" 
case | w |< 1. The same thing happens for the third case, i.e. phantom field. 
We are interested in finding the general solution with the initial conditions 



/(-..«» = /(:). 



t=o 



It is enough to find the fundamental solution of ea. (|52[l A(z,i), satisfying 



A(z,0) = 0, 

since for any set of initial conditions (|54|) we have 



<9A 
~dt 



= S(z) 



tp(z,t) = 



A(z-C,t)/(C) + A(2-C,*)fl(C) 



Indeed 



lim 



<p(z,0) 
d<p(z,t) 

at 



lim 



lim 



d( 



Using a Fourier decomposition for A 



we find the dispersion relation 



j(z). 



A(z,t) 



A(z-C,t)/(C) + A(«-C,*)fl(C) 
^-m 2 A^ f(t) +6(z- 09(0 



/oo 
dke lkz A(k)e 1 ' 
-oo 



U3 



(54) 
(55) 

(56) 

(57) 
(58) 

(59) 
(60) 

(61) 

(62) 



with /i 2 = m 2 / X. Observe that negative kinetic energy, A < 0, behaves as a negative m 2 with respect to propagation. 
Using the initial conditions eas. (|55[) we find 



A(z,t) = 



2tt 



Now we use the integral representation 0) (Integral 6.677.6)] 

1 £ d(;j (^(cty - e) em 2 - z 



sin ct \Jk 2 + /j, 2 

Cy/ k 2 + fl 2 



1 

27 



ai e -ifce 



(63) 

(64) 
(65) 
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where Q[x] is the Heaveside function (Q[x] = 1 for x > and Q[x] — 0, x < 0) and Jo (a;) the Bessel function of zero 
order (notice that the change in integration from eci. (|64[) to ea. l|65l) is taken care by the Heaveside function), to obtain 

A(z,t) = JL y°° dke lkz J J Q (liy/ict)* - e) e-^di (66) 
i roo 

= J dkd(Jo (^V(rf) 2 - e) e[( C t) 2 - ev k{z -^ (67) 

= ij ( M V(ci) 2 -^)e[(ct) 2 -^]. (68) 

The last equation is valid for t 7^ and it gives A(z, i = 0) = 0, as can be seen from ea. (|BT|) . Taking the derivative 
of ea.flfiS'l we have 



rf J x U^{ct) 2 - z< 

A(z,t)= - ^— V , L nUrt^ _ r 

+ 2 Jo (W( ct ) 2 ~ (S(ct -z) + S(ct + z)) (69) 

and for t = we get A(z, 0) = 8(z). Eqs. I|68|l and (|69fl are the main results of this section. They show that a signal 
propagates always with speed smaller or equal than c, since A(z,t) and A(z,t) vanish for \z\ > ct. This ensures that 
a signal generated at t = 0, zo will necessarily arrive at a point \z — Zq\ at a time larger or equal than ct, as can be 
seen from the vanishing of the wave function ip(z, t) for \z — zq\ > ct from eas. (|56fl . (|68|) and l|69|) . 

On the other hand, the result for m 2 < (or A < 0) can be obtained easily from Eq. (|65|l or i|68|l using analytic 
continuation in the parameter /i = ro/vA. In both cases one gets 

Aj(z,t) = lj (l/xlvW 3 ^) 0[(^) 2 - (70) 

where Jo is the modified Bessel function of zero order. The amplitude of ea. (|70|l grows exponentially with time, 
showing explicitly the instability of a field with m 2 < (or A < 0). However, it also shows that the speed of the wave 
packet is smaller or equal to the speed of light, as for the standard scalar field case (m 2 > 0, A > 0). The instability 
growth of the perturbation is ameliorated by the expansion of the universe since instead of growing exponentially 
fast it grows power like. The increase of the amplitude implies that the first order approximation in the perturbed 
differential equation 1)521) will cease to be valid. This implies that we can no longer use ea. (|52|l to study the evolution 
of ip and non-perturbative methods must be used. In this case the universe will be become inhomogeneous. However, 
we would like to emphasis that the sound speed never becomes larger than the speed of light. 

It should be stressed that the above results are general. Indeed, the Principle of Equivalence guarantees that in a 
freely falling reference systems, the equations of motion for the fluid or scalar field will take their Minkowskian form. 

This behavior can be simply illustrated solving ea. (|52ll in 1+1 dimensional space-time. We assume that a rectangular 
semi-infinite wave train of unit amplitude is incident at the spatial point z = and look its subsequent behavior. The 
initial conditions for such a pulse are 

f[z) — sin(kx) * 0(— x) 
dz 

Fig. 2a shows the normal case with fi 2 > and E = u> = 2/i. The front of the train moves towards the right with a 
speed equal to c, but the first crest diminishes and settles slowly into the group velocity. This process is just beginning 
in the figure. The field is stable and there is no increase of the amplitude with time. 

Fig. 2b shows the propagation of the same wave train in the case /i 2 < 0, i.e. when the field is "tachyonic" . In spite 
of that character, the signal velocity is equal to c and the amplitude of the wave vanishes outside the lightcone (i.e. 
for z > t). However, the system is unstable and shows exponential growth of the signal as can be seen from the height 
of the amplitude around z = 0, t = 8. 



V. CONCLUSIONS 



We have shown that scalar fields have naturally regions where \w^ \ > 1 and/or \dp /dp \ > 1 and these properties 
seem to contradict causality. In some of these regions the scalar field can be interpreted as an effective tachyon (i.e. 
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FIG. 2: Propagation of a semi-infinite wave train of frequency ui — 2/i, moving to the right. Fig. a shows the amplitude of the 
pulse as a function of (z, t) in the normal fi 2 > case. Fig.b displays the propagation of the same initial train in the anomalous 
(tachyonic) case fi 2 < 0. In this latter case, the signal shows exponential growth but in both cases the amplitude vanishes 
outside the lightcone, i.e. for z > t. 

m 2 < 0) or phantom (i.e. \<j) 2 < 0) scalar field, i.e. p 2 < in both cases. We studied the evolution of the scalar field 
perturbations for models with \w^ \ > 1 and/or \dp /dp \ > 1 including tachyon, phantom, k-essence and curvature 
models. We have shown that even in these extreme cases the scalar field is consistent with causality and Lorcntz 
invariance provided dispersion is taken into account and the sound speed is never larger than the speed of light. 
Dispersion, however, appears naturally for scalar fields having nontrivial potentials (e.g. mass term) such as those 
used to model dark energy. However, when /i 2 becomes negative the amplitude of the perturbations grows rapidly and 
the first order approximation breaks down. The universe becomes, therefore, inhomogeneous and non-perturbative 
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mathematical methods must be used to study the evolution of the perturbations. 
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